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Abstract.  Regarding  the  present  state  of  knowledge  in  the  field  of  the  turbomachinery 
design,  the  method  for  designing  small  bulb  turbines  with  unequal  specific  work 
distribution  of  the  turbine  runner's  elementary  stages  near  the  hub  is  presented  in  the 
paper.  The  distribution  function  of  specific  work  of  all  the  elementary  stages  is 
obtained,  according  to  which  the  averaged  axisymmetric  flow  surfaces  of  the  turbine 
runner  have  a  negligibly  small  deviation  from  the  cylindrical  flow  surfaces.  The 
specific  work  of  the  near-the-hub  elementary  stages,  in  the  given  distribution  function, 
can  be  reduced  up  to  60%  of  the  required  (design)  specific  work,  still  achieving  nearly 
cylindrical  flow  surfaces. 
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1.  Introduction 

Both  the  hydraulic  turbines  design  and  their  operating  performance  analysis  require 
the  use  of  simpler  methodologies  in  preliminary  design  phases,  especially  when  the 
geometry  of  the  turbine  runner  is  not  completely  defined.  In  the  field  of  the 
turbomachinery  design,  such  attempts  are  continuously  made  thus  providing  us  with 
different  methodologies  and  procedures  for  the  design  process  optimization  [10,  11]. 
With  the  computer  technology  development,  a  significant  breakthrough  of  numerical 
methods  and  numerical  simulation  of  the  fluid  flow  has  been  made,  encouraging  the 
numerical  techniques  incorporation  into  the  design  procedure  and  performance  analysis  of 
the  turbomachinery  [12,  13]. 
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A  number  of  assumptions  and  simplifications,  empirical  equations  as  well  as  the 
designer's  experience  always  accompany  the  designing  process  [1,  2,  5,  7,  9].  The  basic 
assumption  in  the  hydraulic  turbine  designing  procedure  is  axisymmetric  flow  surfaces.  An 
elementary  stage  is  a  flow  space  between  two  elementary  immediate  axisymmetric  flow 
surfaces.  The  intersection  of  stay  vane  (SV),  or  turbine  runner  (TR),  and  axisymmetric  flow 
surfaces,  defines  the  vane  or  blade  profiles  in  the  turbine  elementary  stages. 

Stay  vanes  and  blades  are  formed  by  the  stay  vane  profiles  and  the  runner  blade 
profiles  in  the  turbine's  elementary  stages.  The  vane  and  blade  profiles  in  the  turbine 
elementary  stages  can  be  defined  using  suitable  hydraulic  calculations  for  a  two- 
dimensional  fluid  flow  through  the  profile  cascade  of  axisymmetric  flow  surfaces.  In  the 
cases  when  the  flow  surfaces  are  cylindrical,  the  theory  of  fluid  flow  through  the  straight 
plane  profile  cascade  can  be  used  in  calculations  [6]. 

In  order  to  achieve  cylindrical  or  almost  cylindrical  flow  surfaces  in  the  bulb  turbine 
runners,  these  runners  are  usually  designed  to  obtain  an  equal  specific  work  of  all  the 
elementary  stages.  The  turbine  runners  designed  according  to  such  a  principle  are  made  of 
blade  profiles  with  a  significantly  larger  inclination  angle  near  the  hub  than  near  the  shroud. 

In  order  to  minimize  the  blades'  spatial  curvature  and  to  reduce  the  axial  length  of  the 
runner,  it  is  eligible  to  design  the  runner  with  smaller  specific  work  of  the  elementary 
stage  nearer  the  hub  than  the  blade  periphery  [3,  4].  In  this  way,  the  spatial  curvature  of 
the  stay  vanes  is  also  reduced.  The  possibility  to  design  such  stay  vanes  and  turbine 
runners  is  analyzed  in  the  paper.  Only  small  bulb  turbines,  which  can  be  used  in  small 
hydropower  plants  or  as  models  of  large  water  turbines  (for  the  model  testing  purpose) 
are  considered  in  this  study. 


The  scheme  of  a  bulb  turbine  meridional  cross-section,  with  traces  of  two  elementary 
immediate  axisymmetric  flow  surfaces,  where  one  (below)  is  marked  £1S  Sm,  is  presented 
in  Fig.  1.  A  flow  space  between  two  elementary  immediate  axisymmetric  flow  surfaces 
represents  an  elementary  stage  of  the  runner.  To  determine  the  shape  of  stay  vanes  and 
runner  blades,  it  is  sufficient  (for  small  turbines)  to  define  profiles  in  7  to  12  elementary 
stages,  approximately  evenly  distributed  along  the  blade  height. 


2.  Basic  Formulas 


TR  -  turbine  runner 
SV  -  stay  vane 

Sm  -  flow  surface  of  elementary 


flow  through  the  turbine 
0,  1'-  inlet  and  outlet  control 


cross-section  of  SV 
r(Sm)  -  radius  of  various  control 


cross-section  of  SV 
1,  2-  inlet  and  outlet  control 


cross-sections 


Fig.  1  Scheme  of  a  bulb  turbine  meridional  cross-section 
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The  runner  is  designed  as  a  free-vortex  flow  (to  obtain  zero  circumferential 
component  of  absolute  flow  velocity  at  outlet  (c„2=0));  therefore,  based  on  the  Euler's 
equation,  the  specific  work  of  the  turbine  runner  elementary  stage  is: 

yk(SJ  =  (o-ri(SJ-cul(SJ,  zac„2(SJ=0,  (1) 

where:  cm  cn  cz  -  circumferential,  radial  and  axial  component  of  absolute  flow  velocity 
[m/s]  (absolute  flow  velocity  c=cMu°  +clx°  +c,z° ),  rt  -  radius  in  cross-section  1  [m],  w- 

angular  velocity  of  turbine  runner  [s~  ]. 

According  to  Eq.  (1),  the  required  circumferential  component  of  absolute  flow 
velocity,  in  front  of  the  inlet  of  the  turbine  elementary  stage,  can  be  defined  by  formula: 

cjsm)=  yk(S/:\-  (2) 

Stay  vanes  produce  circumferential  components  of  absolute  velocity  in  front  of  the 
runner  inlet.  Neglecting  the  influence  of  viscous  friction,  the  flow  in  area  between  SV  and 
TR  (vaneless  area)  is  free-vortex  flow: 

r(SJ  ■  cu (SJ  =  const.  =  rv(Sm)- cuV(SJ  =  rx (SJ  ■  cul (SJ , 

therefore,  the  circumferential  components  of  absolute  velocity  of  the  stay  vane  elementary 
stages  can  be  calculated  as: 

cuASJ=^cJSJ  =  ^^-  (3) 

Streamline  inclination  angle  (with  respect  to  circumferential  coordinate  rep  ■  u° ,  where 
u°  is  the  unit  vector  in  the  direction  of  circumferential  velocity  u,  u  =  m-u°,  u  =  r-a>) 
in  the  outlet  of  stay  vane  area  is  calculated  by  the  next  formula: 

a1,(SJ  =  arctgC-^^-,    where  cm=^+c2z  .  (4) 

To  determine  the  inclination  angles  of  vane  profile  chamber  line,  in  the  outlet  of  the 
stay  vane  elementary  stages,  alh(Sm)  =  av(Sm)-Aab  Y(Sm)  ,  it  should  be  assumed  that  the 
value  of  flow  inclination  angle  in  outlet  of  the  stay  vane  is  ( Aabr(Sm)  =  (2-h5)°  ). 

In  front  of  the  inlet  of  the  stay  vane,  there  is  a  zero  circumferential  component  of  the 
absolute  velocity  (cu0  =  0),  thus,  the  inclination  angle  of  the  camber  line  in  the  inlet  of  the 
turbine  stay  vane  is  alM(Sm)  =  90°. 

Design  operating  parameters  of  the  turbine  are:  co,  Q  =  Q+,  YT=  YT+,  where  Q+  and  Yj 
are  the  best  efficiency  (optimal)  operating  parameters  (when  the  turbine  operates  with 
maximum  efficiency,  rj  =  max  77  =  77 +),  Q  [m3/s]  is  turbine  volume  flow  rate  and  YT  [J/kg] 
is  turbine  specific  work  ( YT  =  gHT  (g  =  9.81  m/s2),  HT  -  net  turbine  head  [m]). 

Adopting  the  maximum  value  of  hydraulic  efficiency  ( r/h  =  r/h+),  design  specific  work 
of  the  turbine  runner  is  Yk=  Yk+=  r]h+  Yj  .  When  designing  a  turbine  using  the  model  of 
equal  specific  work  of  turbine  elementary  stages  (yk  (Sm)  =  const.  =  Yk  ),  and  when  flow 
surfaces  in  the  runner  are  cylindrical  (ri(5m)  =  r2(Sm),  czl  =  cz2  =  const.),  the  theory  of  flow 
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in  the  straight  plane  profile  cascades  is  used.  Inclination  angles  of  profile  chord  in 
elementary  stages  of  the  turbine  runner  (/?,  =  pt  (Sm))  are  defined  by  the  lift  force  method, 
where,  using  the  above  mentioned  theory,  the  well  known  formula  is  used  [1,2]: 

^=2^,  (5) 

t  wM 

where:  %y—  profile  lift  coefficient,  I  -  length  of  the  chord  line,  t  -  cascade  spacing, 
Awu  =  wul  -  wu2  (Aw„  =  Ac„,  for  cu2  =  0),  w  =  c-u  is  relative  velocity  of  water  flow, 

w=wau°  +  wrr°  +  wzz°  (  wu  =  cu  -  u,  wr  =  cr ,  wz  =  cz)  and  wm  -  value  of  averaged  velocity 

in  infinity  (w^  =  0.5(w,  +w2)). 

For  elementary  stages  of  the  runner  there  are:  ^v=^v(5m),  t=t(Sm),  l=l(Sm), 
Awu  =  Awu  (Sm)  and  wK  =  wK  (Sm).  Since  coefficient  ^y  depends  on  inclination  angle  of  the 
profile  chord  line,  based  on  formula  (5),  for  calculated  value  ^v,  the  inclination  angle  of 
the  profile  chord  line  can  be  determined  respectively. 

In  the  case  where  flow  surfaces  in  the  runner  are  not  cylindrical  (when  flow 
components  are  not  uniform  along  the  radius),  volume  flow  rate  Q  (Sm)  passing  through 
the  flow  space  between  the  hub  and  flow  surface  Sm,  is  calculated  as  follows: 

1(S,„)  r2{Sm) 

czl(r)rdr  =  2n  j  cz2{r)rdr,  (6) 

where,  czl(r)  and  cz2(r)  are  distribution  functions  of  axial  components  of  absolute  flow 
velocities  in  the  control  cross-sections  in  front  of  and  behind  the  runner.  Functions  cz\(r) 
and  cz2(r)  should  satisfy  the  following  formula  of  volume  flow  through  the  runner: 

re  r 

Q  =  2njczl(r)rdr  =  2njcz2(r)rdr .  (7) 

For  given  function  yk  (r),  the  specific  work  of  turbine  runner  can  be  calculated  using 
the  formula: 

(    r«  ^ 

Y* = i  J y* (r)  • dQ = i  2n\ y* {r)c* {r)rdr  ■  ' 8 '! 


3.  Conditions  for  Obtaining  Cylindrical  Flow  Surfaces  in  the  Control  Cross- 
Sections  in  front  of  and  behind  the  Turbine  Runner 

The  necessary  condition  is  that  the  control  cross-sections  are  placed  in  the  space 
which  is  physically  bounded  by  cylindrical  surface  of  the  hub  (turbine  hub  radius 
r,  =  const.)  and  the  shroud  (turbine  shroud  radius  re  =  const.),  as  shown  in  Fig.  1 . 

The  equation  of  steady  fluid  flow  of  the  inviscid  fluid,  disregarding  gravity  acceleration, 
can  be  written  in  the  form: 
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\c ,  rote  1  =  — -  gradpt ,  (9) 
P 

where  p  is  static  pressure  [Pa]  and  p,  the  total  pressure  ( pt  =  p  +  pc2 12)  [Pa]. 

For  the  fluid  flow  on  cylindrical  surfaces  (cr=  0,  c  =  cuu°  +  c,z° )  in  control  cross-sections 
in  front  of  and  behind  the  turbine  runner  (where  d  I  Sep  =  0  ),  according  to  the  component  of 
vector  Eq.  (9)  in  the  direction  of  coordinate  r,  a  differential  equation  is  obtained: 

c^+c_}Ld(rcJ=ldpL^  (io) 
z  dr     r     dr       p  dr 

which  defines  the  condition  for  cylindrical  flow  surfaces  in  control  cross-sections  1-1  and  2-2. 

In  control  cross-section  2-2  (outlet  of  the  runner)  cll2  =  0  and  p,  =  pt,2  =  const., 
therefore,  according  to  equation  (10)  it  is  c,  =  cz2  =  const. 

Defining  the  hydraulic  efficiency  of  the  runner  elementary  stage  as: 

„  ,r.  ,_p-yk{Sm)  _  p-yk(Sm) 

&P,JSJ  Ph(SJ-pt2(SJ 

and,  assuming  that  hydraulic  efficiency  is  the  same  for  all  elementary  stages 
(%k(Sm)  =  const.  =  t]hk),  the  total  pressure  in  control  cross-section  1-1  (inlet  of  the 
runner)  can  be  expressed  as: 

Ph(r)  =  Ph+^k.  (12) 
r?H.k 

According  to  Eq.  (1): 

rcUl=-yk(r\,    i.e.  c  =  —  yk(r\  (13) 
a>  rco 

where  yk(r\  =  yk{Sm) ,  due  to  rx  =  rt(Sm) . 

Considering  Eqs.  (12)  and  (13),  Eq.  (10)  for  obtaining  cylindrical  surfaces  in  the 
control  cross-sections  1-1  (for  pn(r)  =  pn,  c,  =czl,  cu  =  cu2 )  yields: 


5c„ 


1       1      r  / 

 —  yk(r\ 

V  '\.k    r  ra 


dyk(r\ 


dr 


where:  4^>0. 

For  equal  specific  work  of  the  runner  elementary  stages  (  yk  (r)  =  const  =  Yk  ^  (jue  to 
equation  (14)  it  is  obtained  cz\  =  const.,  and  due  to  Eqs.  (6)  and  (7)  cz\  =  cz2  and  ri(Sm)  = 
r2(Sm),  thus  the  assumption  of  the  cylindrical  flow  surfaces  in  the  runner  is  reasonable  [8]. 

If  dyk  /dr>0,  according  to  equation  (14),  dczi  I  dr  >  0 ,  thus,  due  to  equation  (5),  for 
cz2  =  const.  =  cz ,  can  be  resolved  that  r{{Sm)  >  r2(Sm)  and  flow  surfaces  are  approximately 
conical  in  the  runner. 
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Further  on,  the  relationship  between  the  hydraulic  efficiency  of  turbine  runner 
(f?hk(Sm) )  and  turbine  hydraulic  efficiency  ( %(Sm) )  is  obtained. 

Denoting  ptI  and  ptI1  as  total  pressures  on  inlet  (I)  and  outlet  (II)  of  turbine,  the 
hydraulic  efficiency  of  the  turbine  elementary  stage  is  defined  by  relation: 

^j=^H=  <Ffm\s,-  (15) 

yT(Sm)  P,ASm)-Pt.lliSm) 

Since  Pai^J- ptJ,{Sm)  =       ,  +Ap,  2  +  &ph_n  =  pygI_l  +Pygl_2  +  Pyg2_„  ,  where  ygM 
and  ^si-n  are  specific  losses  of  flow  energy  in  the  flow  passages  from  I  to  1  and  the  flow 
passages  from  2  to  II,  it  is  easy  to  predict  that  there  is  a  relationship  between  t]h(Sm)  and 

Tl*=— I"  or  —,  7'  (16) 

where,  rjh  =%(SJ,  t]hk  =  %.k(SJ,  £,_,  =£,_l(SJ  and  s2_„  =£2_„(SJ  are dimensionless 
losses  of  flow  energy  in  the  passages  from  /  to  1  and  the  flow  passages  from  2  to  // 
(turbine  diffuser), 

,„    ,       yg.I-l(Sm)  ,  .„    .       3^.2-// (^m) 

=  £/-l(SJ  =  777—    and  £2-«  =  W^J  =  777—  •  (16) 

Since  AphJSJ=pyk(SJ  +  pygA_2(SJ  =  p(l  +  el_2(SJ)-yk(Sm),  where  ^_2(SJ  = 
t^s  1-2  (^m  V    (^m ) '  based  on  Eq-  (1 1).  it can  be  written: 

%*(Sm)  =  -  l-7Fr  ■  (17) 

According  to  Eq.  (16),  it  can  be  written: 

%=—  -,  i-e.  s,_l+E1_2+s2_II=—-\.  (17') 

1  +  (^-1+^1-2+  £2-n)  Vh 

Dimensionless  loss  £>_i  contains  flow  energy  losses  of  the  stay  vane.  The  designing  assumption 
is  that  terms  £/_,  and  ex_2  are  larger  than  dimensionless  losses  in  the  diffuser,  £2_n. 


4.  Suggestion  for  Function  yk  (r) 

In  order  to  reduce  specific  work  of  elementary  stages  near  the  hub,  and  to  achieve 
flow  surfaces  inside  the  runner  that  do  not  deviate  much  from  the  cylindrical  surfaces,  the 
distribution  function  of  elementary  stages  specific  work  is  suggested  as  follows: 

yk  (r)  =  A  +  2roB  ■  r  —  B  -r2,  for  rt  <  r  <  r0 
and  yk  (r)  =  yk0  =  const.,       for  rQ  <  r  <  re 

where  is  8yk  (r)/dr  =  0  for  r  =  r0 . 


Design  of  Small  Bulb  Turbines  with  Unequal  Specific  Work  Distribution  of  the  Runner's  Elementary  Stages  79 


Function  yk  (r)  is  defined  according  to  radius  r  in  r[m] 
the  control  cross-section  1-1  (in  front  of  the  runner). 

Function  graph  yk(r),  which  is  defined  by 
equation  (18)  is  presented  in  Fig.  2.  Function  yk(r) 
is  defined  according  to  accepted  parameters  ra  and 
yk.i  =  yk(rd-  To  obtain  flow  surfaces  in  the  turbine 
runner  with  only  a  small  deviation  of  the  cylindrical 
surfaces  the  following  is  recommended: 


ro^-(ri+re)  and 


y,,=  (0.6  +  0.7)7;, 


yk  _yk  Vk(r)[Ikg] 
Fig.  2  Function  yk  (r) 


where  Yk  is  design  specific  work  of  the  runner. 

The  value  of  specific  work  yko  is  an  unknown  value  and  it  can  be  defined  using  a 
requirement  that  specific  work  of  runner  (Yk)  obtained  by  Eq.  (8),  is  equal  to  the  design 
specific  work  of  runner  (Yk  =  Yk). 

Value  yk.o  is  determined,  as  will  be  shown  later,  using  an  iterative  procedure.  The 
coefficients  A  and  B  in  the  first  Eq.  (18)  are  calculated  using  formulae: 


B 


■yu 


(r0-n)2 


and  A=yk( 


(19) 


obtained  from  the  requirement  yk  (/- )  =  yk  i . ,  for  r  =  r,  and  yk  (r0)  =  yk 0 ,  for  r=r0. 

Assuming  that  the  flow  surfaces  in  the  control  cross-section  in  front  of  the  runner  (1- 
1)  are  cylindrical,  based  on  Eqs.  (14)  and  (18),  it  is  obtained: 


dclA 
dr 


a  —  p-r  + 


Y 


S_ 

2 


for  r,  <r<rn 


and  c.j  =  const.  =  czX (r0),    for  rQ  <r  <re 
where  the  coefficients  are: 


(20) 


a  =  4r0B 


'  1 

3B^ 

,  P  =  AB 

'  1  B^ 

+  — 

—  +  — 

.  Y 

CO  j 

4B^fl 

CO 


and  8  ■■ 


4r0AB 
~o7~ 


(21) 


Integral  of  the  first  Eq.  (20),  for  r,  <  r  <  r0,  becomes: 


czi(r)  =  4 Kczi(0)  -a'r0-r)- 


2  r 


o  J 


(22) 


In  the  previously  shown  function  czl(r),  for  r,  <  r  <  r0 ,  there  is  an  unknown  velocity 
value  c-i(r0)  (axial  component  of  flow  velocity  cz\  =  czl(r0),  for  rQ<  r<  re).  This  velocity 
can  be  defined  according  to  the  requirement  that  the  flow  rate  calculated  by  Eq.  (7)  is 
equal  to  the  design  turbine  flow  rate  (Q  =  Q+).  Since  coefficients  A  and  B,  and  also 
coefficients  a,  /?,  y  and  S,  depend  on  unknown  variable  y^o,  it  can  be  concluded  that 
function  c-i(r0)  is  determined  using  an  iterative  procedure. 
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Since  czl(r)  =  czl(r0),  for  r0<  r  <  re,  and  for  r,<  r  <  r0  czl(r)  is  defined  by  Eq.  (22),  Eqs. 
(7)  and  (8),  used  for  determination  of  values  yk0  and  czl(r0),  can  be  derived  into  forms: 

Q  =  2njczl(r)rdr  +  nczl(r0)(r>-r>)  (23) 

Y*=^  M  %  ^  ^rdr + ^  ^^.oc,i  WO;2  -  '•o2))  •  (24) 


5.  Determination  of  yk0  and  czl(r0) 

Values  0  and  czl(r0)  are  determined  by  using  an  iterative  procedure.  Perhaps,  it  is 
better  to  say  that  this  is  a  double  iterative  procedure,  since  in  each  step  of  determining 
yk  0 ,  czl(r0)  should  be  determined  as  well  by  another  iterative  procedure  [3]. 

In  order  to  reduce  the  number  of  iterative  steps  for  determining  yk0,  this  value  is 
obtained  in  the  first  iterative  step,  using  a  formula: 

ya=y:*:y"'  (25) 

l-k 

where  k  =  k  (r0,  rh  re)  is: 

k  =  j      -n)  +  0-5^  ~  r4  )  - 1 .333r0  (r03  -  r3 ) 

(';2-^2)('-o-';)2 

Formula  (25)  is  derived  using  the  assumption  that  the  axial  components  of  flow  velocities  in 
the  control  cross-section  in  front  of  the  runner  change  negligibly  (for  czl  =  const.  =  cz ). 

Using  value  y{kl\ ,  which  was  determined  by  formula  (25),  the  procedure  of  determining 
yto  is  completed  in  maximum  of  three  iterative  steps.  The  reason  for  this,  as  calculations 
showed,  is  that  czl(r)  values  change  less  than  specific  work  yk  (r),  for  r,  <  r  <  r0.  In 
iterative  steps  of  determining  yk0  correction  of  this  value  is  performed,  requiring  that 
calculated  specific  work  (Yk)  using  formula  (24)  slightly  differs  from  the  design  specific 
work  of  the  turbine  runner  Yk  (\Yk-  Y*\lY*  <  0.005 ). 

In  each  iterative  procedure  of  determining  yL0  the  value  czi(r0)  is  also  obtained  by  an 
iterative  procedure.  In  the  first  iterative  step  it  can  be  taken  that  czl(r0)  =  1.2c. ,  where: 

ct=Q/[x(r?-r?)],  for  Q  =  Q+.  (26) 

In  iterative  steps  of  determining  czl(r0)  a  correction  of  this  value  is  performed,  and 
according  to  formula  (23),  calculated  flow  rate  Q  slightly  differs  from  design  flow  rate  Q+ 

(\Q-Q+\iQ+<  0.005). 

In  the  program  for  determining  y^.o  and  czl(r0)  [4],  in  the  iterative  procedure  of  solving 
the  problem,  czl(r0)  changes  with  the  iterative  step  +0.005  •  c  ,  and  yk0  changes  with  the 
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step  +0.005  •  Yk+  .  The  integrals  in  formulae  (23)  and  (24)  are  calculated  using  a  trapezoid 

rule,  whereas  in  the  integration  area  from  r,  to  r0  calculation  points  are  distributed  for  the 
iterative  step  Ar  =  0.0025m  (2.5mm).  Besides  printing  the  calculated  values  y^.o  and 
czi(roX  program  [4]  is  also  developed  to  print  values  yk{r),  czl(r)  and  Q(r),  for  r  e  [r;,rj, 
for  each  iterative  step  Ar  =  0.0025m,  where  r  -  radius  in  control  cross-section  in  front  of 
the  turbine  runner  and  Q(r)  is  the  volume  flow  rate  under  the  cylindrical  flow  surface  of 
radius  r. 

In  the  control  cross-sections  before  (1-1)  and  behind  (2-2)  turbine  runner,  the  flow  surfaces 
are  cylindrical,  where  cz2  =  const.  =  cz .  According  to  Q(r)  data,  where  r  =  n(5m),  the  functional 
relationship  of  cylinder  radius  in  the  same  axisymmetric  flow  surfaces  can  be  obtained. 
Regarding  equation  (6),  for  Q(Sm)  =  Q(r),  r2(Sm)  =  r2(f)  and  cz2  =  const.  =  cz  it  follows: 


where  r  =  rx (Sm )  and  r2  =  r2(SJ. 

For  r2(r)  <  1,  as  obtained  for  dyk  /  dr  >  0 ,  the  flow  surfaces  in  the  turbine  runner  are 
approximately  conical  and  n(5m)  =  r2(Sm). 

To  determine  the  real  flow  surfaces'  deviation  from  the  cylindrical  shape,  the 
dimensionless  function  is  1/2(r)  =  r I r2(r)  ■  In  cases  where  rl/2(r)  <  1.03,  it  can  be  said  that 
the  real  flow  surface  in  the  turbine  runner  negligibly  vary  from  the  cylindrical  surface. 

The  same  calculation  procedure  has  been  used  for  determining  values  yk0  and  cz2(r0) 
and  coefficients  A,  B,  a,  p,  ywd  <51ike  in  the  computer  program  given  in  Ref.  [4].  Even 
though  the  program  is  originally  created  for  axial  flow  fans,  the  same  code  can  be  used  for 
bulb  turbine  as  well,  replacing  value  t]h  with  l/r|h  and  value  cz2(r0)  with  cd(r0). 

Diagrams  of  yk(r),  czl(r)  and  r1/2(r)  are  given  in  Fig.  3,  for  ykj  =  0.7Tj.+  and  different 
values  of  r0  (r0=  180,  210,  240,  280  and  310  mm)  of  one  small  bulb  turbine,  which  is  built  in 
the  small  hydro  power  plant  "Grcki  mlin"  near  Prokuplje.  The  geometrical  and  design 
operating  parameters  of  the  turbine  are:  r,=  142  mm,  re=355  mm,  Q+=  1.4m3/s, 
Yk+=  20.36  J/kg  and  co  =  41.9  r/s  (n  =  400  min1).  For  adopted  r\  =  0.75  (rjm=  0,94,  t]Q=  0.96 
and  £i_2=  +  £2.n)  it  follows  that  rjhk=  0.91.  As  can  be  concluded  from  Fig.  3,  for 
yu=  0.71V"  =  14.2  J/kg,  r1/2(r)  <  1.03  is  obtained  (when  flow  surfaces  in  the  runner  can  be 
considered  cylindrical),  for  r0<210mm.  It  can  be  shown  easily  that  for  values 
ykJ  =  0.6Yk+=  12.2  J/kg,  rV2(r)  <  1.03  is  obtained  for  r0  <  190  mm.  For  the  observed  turbine 
it  is  n  I  re  =  0.40. 

The  bulb  turbine  in  the  small  hydro  power  plant  "Grcki  mlin"  on  the  river  Toplica, 
near  Prokuplje,  is  a  small  turbine  (PT=  26  kW);  therefore,  the  blades  are  shaped  using 
panel-like  profiles  of  the  runner  elementary  stages.  The  runner  consists  of  4  blades,  and 
due  to  the  technical  reasons  the  blade  profile  thickness  is:  <$■=  16  mm  near  the  hub  and 
Se=  10  mm  on  the  blade  periphery.  Due  to  structural  constraints,  the  profile  lengths  are 
/,  =  296  mm  (r,  =  142  mm)  and  le  =  525  mm  (re  =  355mm).  The  requirement  is  to  maintain 
geometric  parameters  r,  =  142  mm  and  rf=  355mm. 


(27) 
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The  calculation  based  on  the  model  of  unequal  specific  work  of  the  elementary  stages, 
that  are  changed  by  the  rule  given  in  equation  (18),  for  yki=0.6Yk+  =  12.2  J/kg  and 
rQ  <  180  mm  (when  cz_,(r0)  =  4.27  m/s  ,  yk0  =  20.68  J/kg,  A  =  -169.5,  B  =  5871,  a=  47065, 
p=  104426,  r=  "7356,  8=  -408.2  and  i[/2(r)  <  1.023  )  the  following  is  obtained:  /?,,=  31° 
and  pte=  15.8°.  The  inclination  angle  of  the  blade  profile  is  /?,,  -  pte  =  15.2°,  which  is  6.2° 
(29%)  smaller  than  the  blade  profile  designed  based  on  the  model  of  the  equal  specific 
work  of  the  turbine  elementary  stages. 

The  values  of  iterative  steps  and  defined  relative  errors  of  specific  work  and  volume 
flow  rate  can  be  selected  in  the  program  for  determination  of  yk0  and  czl(r0)  [4].  However, 
the  values  applied  in  the  above  mentioned  example  of  the  bulb  turbine  give  the  results 
that  are  accurate  enough  for  the  technical  practice. 


r=n  [mm] 
re=355 


r=n  (Sm)  [mm] 


r,(S.) 


10   12   14   16    18   20  22 


»yt  [J/kg] 
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Fig.  3  Diagrams  of  yk(r),  czl(r)  and  rl/2(r) 


Design  of  Small  Bulb  Turbines  with  Unequal  Specific  Work  Distribution  of  the  Runner's  Elementary  Stages  83 

Conclusion 

The  distribution  function  of  specific  work  in  turbine  elementary  stages  {ydj))  are 
proposed  in  the  paper.  In  addition,  function  ri/2(r)  is  defined,  according  to  which  the 
flow  surfaces  deviation  compared  to  the  cylindrical  shapes  can  be  determined. 

Using  the  proposed  functional  distribution  of  specific  work  of  the  turbine  elementary 
stages,  with  the  fulfilled  condition  that  the  flow  surfaces  negligibly  deviate  from  the 
cylindrical  surfaces,  the  designed  bulb  turbine  blades  are  less  twisted.  This  design  method 
is  applicable  to  stay  vanes,  obtaining  less  twisted  stay  vanes  as  well. 

Acknowledgement:  This  paper  is  result  of  technological  project  No.  TR33040,  "Revitalization  of 
existing  and  designing  new  micro  and  mini  hydropower  plants  (from  100  to  1000  kW)  in  the 
territory  of  South  and  Southeast  Serbia  "  supported  by  Ministry  of  Science  and  Technological 
Development  of  the  Republic  of  Serbia. 


References 

1.  Babic,  M.,  Stojkovic,  S.,  1990,  Basics  of  Turbomachinery:  operating  principles  and  mathematical 
modeling,  University  of  Kragujevac,  Naucna  knjiga,  Beograd,  Serbia  (in  Serbian). 

2.  Benisek,  M.,  1998,  Hydraulic  Turbines,  Faculty  of  Mechanical  Engineering,  Belgrade  (in  Serbian). 

3.  Bogdanovic,  B.,  Bogdanovic-Jovanovic,  J.,  Spasic,  Z.,  2010,  Designing  of  low  pressure  axial  flow  fans  with 
different  specific  work  of  elementary  stages,  The  International  Conference  Proceedings:  Mechanical 
Engineering  in  XXI  Century,  pp.  99-102,  Nis. 

4.  Bogdanovic,  B.,  Bogdanovic-Jovanovic,  J.,  Todorovic,  N.,  2011,  Program  for  determination  of  unequal 
specific  work  distribution  of  elementary  stages  in  the  low-pressure  axial  flow  fan  design  procedure,  Facta 
Universitatis,  Series  Mechanical  Engineering,  9(2)  pp. 149-160. 

5.  Obradovic,  N.,  1974,  Turbocompressors,  Tehnicka  knjiga,  Beograd  (in  Serbian). 

6.  Bogdanovic,  B.,  Bogdanovic-Jovanovic,  J.,  Spasic,  Z.,  Milanovic,  S.,  2009,  Reversible  axial  fan  with  blades 
created  of  slightly  distorted  panel  profiles, ,  Facta  Universitatis  Series  Mechanical  Engineering,  7(1)  pp.23-36 

7.  Ristic  B.,  Milenkovic,  1996,  Small  hydropower  plants,  Naucna  knjiga,  Beograd  (in  Serbian). 

8.  Bogdanovic-Jovanovic,  J.,  Bogdanovic,  B.,  Milenkovic,  D.,  2012,  Determination  of  averaged  axisymmetric 
flow  surfaces  according  to  results  obtained  by  numerical  simulation  of  flow  in  turbomachinery, 
Thermal  Science,  16(2)  pp.  647-662. 

9.  Jovicic,  N.,  Babic,  M.,  Jovicic,  G.,  Gordic,  D.,  2005,  Performance  Prediction  of  Hydraulic  Turbomachinery, 
Facta  Universitatis,  series  Mechanical  Engineering,  3(1)  pp.  41-57. 

10.  Albuquerque,  R.B.F.,  Manzanares-Filho,  N.,  Oliveira,  W.,  2007,  Conceptual  optimization  of  axial-flow 
hydraulic  turbines  with  non-free  vortex  design,  IMechE  221,  Part  A:  J.  Power  and  Energy,  pp. 713-725. 

11.  Hofler,  E.,  Gale,  J.,  Bergant,  A.,  201 1,  Hydraulic  design  and  analysis  of  the  saxo-type  vertical  axial  turbine, 
Transactions  of  the  Canadian  Society  for  Mechanical  Engineering,  35(1),  pp.1 19-143. 

12.  Da  Cruz,  A.G.B.,  Luiz,  A.  Mesquita,  A.L.A.,  Blanco,  C.J.C,  2008,  Minimum  Pressure  Coefficient  Criterion 
Applied  in  Axial-Flow  Hydraulic  Turbines,  J.of  the  Braz.Soc.  of  Mech.  Sci.  &  Eng.,  30(1),  pp. 30-38. 

13.  Ferro,  L.M.C.,  Gato,  L.M.C.,  Falcao  A.F.O.,  2011,  Design  of  the  rotor  blades  of  a  mini  hydraulic  bulb- 
turbine,  Renewable  Energy  36,  pp.2395-2403. 


84 


J.  BOGDANOVIC-JOVANOVIC,  B.  BOGDANOVIC,  I.  BOZIC 


PROJEKTOVANJE  MALIH  CEVNIH  TURBINA  SA  RAZLICITIM 
JEDINICNIM  RADOVIMA  ELEMENTARNIH  STUPNJEVA 

OBRTNOG  KOLA 

Uzimajuci  u  obzir  dosadasnja  saznanja  iz  oblasti  projetovanja  turbomasina,  u  ovom  radu  je 
prikazan  metod  projektovanja  malih  cevnih  turbina  sa  razlicitim  jedinicnim  radovima  elementarnih 
stupnjeva  turbinskog  kola  u  okolini  glavcine.  Data  je  funkcija  raspodele  jedinicnih  radova 
elementarnih  stupnjeva  turbinskog  kola,  pri  kojoj  osnosimetricne  strujne  povrsine  u  turbinskom  kolu 
zanemarljivo  malo  odstupaju  od  cilindricnih  strujnih  povrsina.  U  datoj  funkciji  raspodele,  jedinicni 
rod  elementarnog  stupnja  turbinskog  kola  moze  se,  uz  glavcinu,  smanjiti  i  na  60%  proracunskog 
jedinicnog  rada  turbinskog  kola,  a  da  strujne  povrsine  u  obtnom  kolu  budu  prilizno  cilindricne. 

Kljucne  reci:  cevna  turbina,  elementarni  stupanjevi,  jedinicni  rad,  osnosimetricne  strujne  povrsine. 


